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UCUMITV  CLASSIFICATION  OP  THIS  PASinNb'a 


STARK  PROFILE  CALCULATIONS  FOR  LTMAM  SERIES  LIMES 


OF  ONE-ELECTRON  IONS  IN  DENSE  PLASMAS 

I.  INTRODUCTION 

Interest  In  Stark  broadened  spectral  lines  emitted  or  absorbed  by 
dense  high  temperature  plasmas  arises  mainly  from  three  classes  of 
problems.  In  conjunction  with  a reliable  theory  of  the  broadening, 
measured  profiles  can  first  be  used  to  determine  one  of  the  most  Impor- 
tant parameters  of  the  plasma,  namely  Its  density.  The  theory  applied 
In  the  present  work  actually  yields  profiles  of  absorption  and  emission 
coefficients  (or  cross  sections)  for  plasmas  In  which  Stark  broadening 
due  to  thermal  fluctuations  of  the  electric  mlcrofleld  dominates  all 
other  line  broadening  mechanisms,  including  Stark  effects  arising  from 
plasma  waves  and  microinstabilities.  The  corresponding  restrictions 
must  be  kept  in  mind  when  applying  this  method  of  density  diagnostics, 
e.g.,  to  plasmas  used  in  laser  or  particle-beam  fusion  research 
(Imploding  pellets). 

A second  major  class  Is  the  calculation  of  radiative  transfer  in 
stellar  Interiors  and,  again,  imploding  pellet  plasmas.  For  such  appli- 
cations, one  is  mainly  interested  in  the  behavior  of  absorption  coeffi- 
cients at  large  frequency  separations,  compared  to  the  half  width,  from 
line  center.  Also,  since  Stark  profiles  follow  approximately  an  inverse 
2.0  - 2.5  power  law  as  function  of  this  separation.  Stark  broadening  may 
be  dominant  in  this  context  even  in  circumstances  where  the  corresponding 
half  width  is  smaller  than  the  Doppler  width.  This  can  happen  because 

thermal  Doppler  profiles  are  Gaussian,  l.e.,  decay  exponentially  on  the 
wings  of  the  lines. 


Note:  Manuscript  submitted  July  17,  1978. 
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The  third  class  of  problems  is  concerned  with  turbulent  plasmas 


in  which  fields  from  strongly  excited  collective  modes  may  have  a sub 


stantlal  influence  on  spectral  line  shapes.  Appropriate  measurements 


and  analysis  enable  one  to  determine  electrical  energy  densities  and 


dominant  frequencies  in  such  cases.  We  shall  assume,  however,  that  wave' 


produced  fields  are  of  minor  importance  in  comparison  with  particle' 


produced  fields,  and  merely  refer  the  reader  to  an  Introduction  to  this 


rapidly  developing  method  of  plasma  diagnostics 


Particle-produced  fields  are,  on  the  average,  a factor  ~3  larger 


than  the  Holtamark  normal  fleldstrength 


2.603peN' 


This  is  very  nearly  equal  to  the  field  produced  by  an  ion  of  charge  p 

which  is  at  a mean  ion-ion  distance  r « (^ffN  /3)  from  the  perturbed 

P P 

ion.  Since  in  a neutral  plasma  containing,  e.g.,  only  a single  ion 

species,  the  electron  density  is  related  to  the  ion  density  by  N = pN  , 

P 

the  electronic  Holtsmark  field  , 


is  related  to  that  produced  by  ions  by 


Except  for  very  highly  charged  ions,  typical  fields  produced  by  the  two 
kinds  of  charged  particles  are  therefore  of  the  same  order,  and  we  may 
u*®»  •»y»  5*0  f°r  estimating  Stark  effects  in  thermal  plasmas. 

To  the  extent  that  the  actual  time  dependence  of  the  perturbing 


fields  can  be  neglected  (quaaistatic  approximation,  see  the  following 
section),  linear  Stark  shifts  are  therefore  of  the  order^  (in  angular 
frequencies) 


2 2 
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where  a^  ■ h /me  is  the  Bohr  radius  and  n the  principal  quantum  number 
of  the  upper  state  of  the  radiating  ion  with  nuclear  charge  z.  For  the 

Stark  effect  to  be  linear  rather  than  quadratic.  Am  e must  be  much  larger 

s 

than  the  fine  structure  spread^  of  the  upper  states,  l.e., 

2 4 

Am  g » Am  F « 3£_  (n-1)  (5) 


(Here  a w 1/137  1*  the  fine- structure  constant.)  He  may  define  a 
critical  electron  density  Np  by  the  corresponding  equality.  This  den- 
sity is 


w or  m 


\15/2 


A3/2 


(6) 


and  must  be  well  exceeded  for  our  calculations  to  be  valid. 

Another  characteristic  density,  N^,  is  obtained  by  comparing  Am  g 

with  the  thermal  Doppler  width,  which  la  for  typical  temperatures 

12  2 

kT  m •g  me  (az)  and  ion  masses  m^  « 2lEmH  estimated  by 


(7) 


For  density  measurements  based  on  our  calculations,  a second  lower  limit 
is  therefore  given  by 


5 21A 


3 18  «*3 

again  expressed  in  terms  of  (a/a^)^  « 3 x 10  cm  . We  also  note 

2 q/k 

N_/N  « (n  /z)  , which  indicates  that  fine  structure  is  more  important 

1)  F 

2 

than  Doppler  broadening  for  n < z.  (Natural  broadening  need  not  be 
considered  because  it  is  always  less  than  the  fine-structure  splitting.) 

High  density  limits  for  the  validity  or  usefulness  of  the  present 
calculations  might  be  expected  from  a number  of  sources.  First  of  all, 
for  the  lines  of  different  principal  quantum  numbers  not  to  overlap 
and  the  linear  Stark  effect  approximation  to  remain  valid.  Stark  shifts 
should  certainly  be  smaller  than  separations  between  the  upper  level  in 
question  and  the  next  level, 

mc^  . . 2 T 1 1 1 a 2.  .2  -3 

•Wl-"n  - ST  (OT)  [7  - ^2 1 “ * c (c“)  - • <»> 

Near-equality  of  this  separation  with  Atu  g from  Eq.  (4)  corresponds  to 


the  Inglis-Teller  limit,  in  our  case 


-2  -A  t/'  £ 

Nn«2xlOia03  * 1.5  * 10 

n 


As  can  be  seen  from  NIT/ND  M 10^z"^\T'^2  and  NIT/NF  * 10^z-^,  there 
is  a considerable  range  (by  factors  10^  to  10^)  between  the  lowest  den- 
sities of  Interest  and  the  Inglis-Teller  limit  for  the  lines  and  ions 
considered  In  this  paper.  However,  it  is  not  obvious  that  this  entire 
range  can  be  explored  with  the  usual  approximations  of  Stark  broadening 
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theory,  which  will  be  discussed  next.  Limitations  Imposed  by  the 
various  approximations  are  further  discussed  In  Section  III  and  in 
Section  IV,  where  we  present  quantum-theoretical  results  for  the 
broadening  by  electrons.  Stark  profiles  and  line  widths  are  given  and 
discussed  in  Section  V. 

II.  THEORY 

As  in  many  previous  calculations  of  Stark  broadening, ^ we  use  the 
following  expression  for  the  spectral  line  shape, 

00 

Ui»)  = - - 8.  Tr  f dF  W(F)  lCF+»f ')  . (11) 

’ Jo 

Here  &u>  is  the  frequency  separation  from  the  unperturbed  line.  Re  in- 
dicates the  real  part  and  Tr  the  trace  over  unperturbed  states  of  the 
radiating  ion  for  which  the  principal  quantum  number  equals  that  of  the 
upper  state  of  the  line  in  question.  The  integral  is  over  the  field 
strength  F produced  by  the  perturbing  ions,  th  W(F)  being  the  corre- 
sponding distribution  function.'’  Hie  dipole  operator  D generates  appro- 
priate products  of  matrix'  elements  between  states  of  principal  quantum 
number  n and  the  ground  state,  and  C is  an  operator  whose  matrix  elements 
(in  terms  of  parabolic  wave  functions)  are  the  linear  Stark  coefficients 
for  the  components  of  the  line.  Finally,  0 is  an  operator  describing 
the  effects  of  electron  collisions  on  the  line  shape. 

To  the  leading  orders  in  both  the  perturbation  (Dyson  series)  ex- 
pansion for  the  (perturbing)  electron- (radiating)  ion  system  and  the 
multipole  expansion  for  the  interaction  Hamiltonian,  this  collision 
operator  is1*^ 

-5- 


f 


(12) 


where  v is  the  velocity  of  the  perturbing  electrons,  and  R is  the 

(radiator)  electron  position  operator  (in  atomic  units)  operating  in  the 

subspace  of  principal  quantum  number  n.  The  quantities  0 and  n 

Kmax  Main 

correspond  to  the  limits  of  the  integral  over  impact  parameters,  which 
arises  in  the  straight  classical  path  approximation  (for  the  perturbing 
electrons)  used  in  Ref.  6.  (The  Impact  parameters  p were  of  course 
assumed  to  be  larger  than  the  radial  coordinate  of  the  radiating 
electron.) 

At  small  impact  parameters,  curvature  of  the  classical  path  due 
to  the  strong  Coulomb  interaction  with  the  radiating  ion  may  be  important. 
Because  the  dipole-monopole  Interaction  can  now  change  sign,  a hyperbolic 

7 

classical  path  calculation'  actually  gives  a convergent  result  fra-  which 

an  equivalent  (Coulomb)  cutoff  follows  as  pc  « (z-l)e2/mv2.  However, 

since  this  is  smaller  than  the  excited  state  Bohr  radius,  0 « 2a,Vz, 

pn  n 0 

^ _ 2 2 2 

by  a factor  ~ (ze  /n  ftv)  <1,  other  than  dipole  terms  in  the  multipole 
expansion  (and  penetrating  orbits)  will  be  very  Important  for  p w P(}. 
Higher  order  terms  in  the  perturbation  expansion  must  be  considered  as 

well,  giving  rise  to  an  equivalent  (strong  collision)  cutoff 

2 

pfl  w n fi/zav  in  the  straight  classical  path  calculations,  to  be  supple- 
mented by  a quantum-mechanical  cutoff  « fi/tnv  corresponding  to  the 
DeBroglie  wavelength  of  the  perturbing  electron.  Because  of 

p J pn  = a /Av  < 1 for  our  high  temperature  plasmas  and,  even, 

2 2 

Pq/p„  * *•  /n  fiv  < 1 we  conclude  that  pR  is  always  the  largest  minimum 
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Impact  parameter  for  typical  electron  velocities  in  plasmas  emitting  the 
lines  under  consideration.  We  therefore  provisionally  choose 


Pmin  ~~  z a0 


Corrections  and  errors  associated  with  this  choice  will  be  dis- 
cussed in  the  following  section,  and  we  now  proceed  with  the  determination 

1 8 

of  the  mn-r-imiiTn  impact  parameter,  n . Three  physical  effects  ’ were 
ignored  in  arriving  at  the  divergent  expression  for  the  collision 
operator:  correlations  between  perturbing  electrons,  splitting  of  the 
levels  by  quasistatic  fields,  and  the  finite  duration  of  the  collisions. 
The  equivalent  cutoff  accounting  for  electron-electron  correlations  is 
close  to  the  (electron)  Debye  radius. 


/ kT  \ 1/2  v_ 

\4flNe2  / U)P  ' 


2 1/2 

where  <Dp  = (4flNe  /m)  is  the  electron  plasma  frequency.  The  quasi- 
static  splitting  may  be  allowed  for  by  p « v/a<dc  with  A<b_  estimated 
by  Eq.  ( k ),  while  the  finite  duration  of  the  collisions  requires 
p <v/|&u)|.  The  combined  effects  are  accounted  for  by  choosing 


2 2 2-1/2 
Pmax  = V(“P  + + A®  > 


and  the  corresponding  theoretical  uncertainties  will  again  be  discussed 
in  the  following  section.  (Since  deviations  from  straight  paths  are  nbt 
Important  near  p^,  we  may  neglect  them  for  pM^.) 

From  Eqs.  (12),  (13)  and  (15)  the  Maxwell  average  of  the  collision 
operator  is  as  in  Eq.  (29)  of  ref.  6 


1 


(16) 


with 


y = 


2 / ov  2 A 2 . 2 

I min  /ftnf\  “V  4 *“s  + *» 

2kT  * \2z  / ^kT 


(17) 


determined  from  the  requirement  pmOT  > Pm^n*  C^e  use  ^ = e ^ao*  *,e,» 
the  Ionization  energy  of  hydrogen,  to  combine  the  various  atomic 
constants.)  The  fraction  of  electrons  excluded  by  this  requirement  is 
estimated^  by  j ir-1^2y^2  and  must  be  small  for  our  calculations  of 
the  broadening  by  electrons  to  be  reasonably  accurate. 

The  most  important  quantity  that  remains  to  be  determined  Is  the 
ion  field  strength  distribution  function  W(F)  in  Eq.  (11)  for  the  line 
shape.  We  have  used  distribution  functions  calculated  by  Tighe  and 
Hooper,^  which  depend  on  electron  and  ion  densities  and  on  the  ion 
charges.  It  is  convenient  to  use  reduced  field  strengths  defined  by 


F/F0 

with  Fq  according  to  Eq.  (2).  The  0 distribution,  besides  on  ionic 
charges  and  ion  density  ratios,  depends  on  the  dimensionless  parameter 


This  parameter  is  primarily  a measure  of  how  strongly  coupled  singly 

2 2 

charged  particles  are.  (Note  that  a /3  = e /rQkT  is  .ue  ratio  of 
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characteristic  potential  and  kinetic  energies.)  The  parameter  a is 


III.  ERRORS  AND  RESTRICTIONS 

The  two  parameters  in  our  problem  that  should  be  small  are  the 
quantities  y for  the  electron  broadening  and,  at  least,  a1/,2/p  for 

the  ion  broadening.  They  are  both  increasing  functions  of  density  and 

decreasing  functions  of  temperature.  For  any  given  ion,  the  temperature 

1 2 

range  is  relatively  small,  and  we  shall  simply  use  kT  = jj-  z EH  in  this 
section.  Also,  the  Inglls-Teller  estimate  should  give  a reasonable  upper 
density  limit  for  the  lines  to  remain  discrete,  provided  0.1  ^ 

is  Indeed  small.  To  verify  this  condition,  we  substitute  NIT  from 
Eq.  (10)  into  Eq.  (19)  and  obtain 


also  a measure  of  the  relative  contributions  of  wave-  and  partlcle- 

1 1/2,— — 

produced  fields,  their  ratio  being  estimated  by  0.17e  / p , if  p 

is  the  mean  charge  of  perturbing  ions.  Since  only  static  correlations 

are  allowed  for  in  the  calculated  field  strength  distribution  functions, 

we  must  certainly  restrict  our  calculations  to  a^2/p*^  < 1.  However, 

— -1/2 

since  for  fixed  electron  density  the  ion  Debye  radius  scales  as  (p) 

and  since  the  mean  ion-ion  separation  is  p^^rQ>  the  generalized  para- 

_5/6 

meter  corresponding  to  a is  a factor  p larger  than  the  value  given 
by  Eq.  (19) . For  the  modified  cluster  expansion  used  in  the  calculation 

of  the  distribution  function  to  converge  rapidly  one  would,  a priori, 

_5/6- 

therefore  expect  the  more  severe  constraint  ap  <1. 


. < 5 CnV1'4  < ,-1A 
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characteristic  potential  and  kinetic  energies.)  The  parameter  a is 


also  a measure  of  the  relative  contributions  of  wave-  and  partlcle- 

1 1/2  —^-5  — 

produced  fields,  their  ratio  being  estimated  by  0.17a  / p , if  p 

is  the  mean  charge  of  perturbing  ions.  Since  only  static  correlations 

are  allowed  for  in  the  calculated  field  strength  distribution  functions, 

we  must  certainly  restrict  our  calculations  to  a^^/p1^  < 1.  However, 

— -1/2 

since  for  fixed  electron  density  the  ion  Debye  radius  scales  as  (p) 

and  since  the  mean  ion-ion  separation  is  p*^rQ,  the  generalized  para- 

5/6 

meter  corresponding  to  a is  a factor  p larger  than  the  value  given 
by  Eq.  (19).  For  the  modified  cluster  expansion  used  in  the  calculation 

of  the  distribution  function  to  converge  rapidly  one  would,  a priori, 

_5 16- 

therefore  expect  the  more  severe  constraint  ap  <1. 


III.  ERRORS  AND  RESTRICTIONS 

The  two  parameters  in  our  problem  that  should  be  small  are  the 
quantities  y for  the  electron  broadening  and,  at  least,  a^^/p  for 

the  ion  broadening.  They  are  both  increasing  functions  of  density  and 
decreasing  functions  of  temperature.  For  any  given  ion,  the  temperature 
range  is  relatively  small,  and  we  shall  simply  use  kT  = ^ z^E^  in  this 
section.  Also,  the  Inglls-Teller  estimate  should  give  a reasonable  upper 
density  limit  for  the  lines  to  remain  discrete,  provided  0.17a^^/p  ^ 
is  indeed  small.  To  verify  this  condition,  we  substitute  from 
Eq.  (10)  into  Eq.  (19)  and  obtain 

. , . 5N-l/4  . -1/4 

a < 3 (zn')  < z 
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The  omitted  wave  contribution  to  the  mlcrofleld  are  therefore  for  p m z 

of  order  0.2  tfU2h  or  less,  and  we  may  conjecture  that  corresponding 

theoretical  uncertainties  in  the  ion  broadening  calculations  are  < 10^ 

for  the  Lyman-a  lines  (and  the  highest  densities)  and  < 5$  for  the  other 

lines.  Should  the  wave  fields  exceed  their  appropriate  thermal  levels 

or  act  mostly  on  a small  portion  of  the  line  profile,  e.g.,  near  the 

center  of  Lyman  a10*11  (see  also  Section  V),  their  omission  here  would 

cause  proportionally  larger  errors.  Equally  serious  appears,  especially 

for  the  o-lines,  the  possible  violation  of  the  constraint  ap  < 1. 

12 

However,  comparison  with  Monte  Carlo  calculations  suggests  that  the 

calculated  distribution  functions  remain  accurate  to  better  than  10£  up 

5/6 

to  az  m2.  The  corresponding  limit  on  the  density  is  often  more 

restrictive  than  that  imposed  by  approximations  in  the  electron  broadening. 

Other  errors  in  the  ion  broadening  calculations  are  connected  with 

the  use  of  the  linear  Stark  effect,  dipole  interaction  approximations. 

While  quadratic  Stark  effects  should  be  small,  except  near  the  Inglls- 

Teller  limit,  quadrupole  interactions  are  Important  already  at  lower 

15 

densities,  as  emphasized  by  Demurs  and  Sholln.  These  interactions 
cause  asymmetries  in  the  line  shapes  of  order 


if  the  Inglis-Tellar  limit  is  used  for  an  upper  bound.  Especially  for 
the  Lyman-a  lines,  one  might  therefore  expect  substantial  asyonetries 
from  this  source.  However,  detailed  calculations1^  of  the  ion-produced 
asymmetry  yield  additional  factors  < 0.5  so  that  omission  of  the  higher 
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order  effects  discussed  In  this  paragraph  should  generally  be  less 
serious  than  other  uncertainties  in  our  calculations. 

Before  proceeding  to  the  electron  broadening,  we  wish  to  suggest 
that  errors  introduced  by  the  quasistatic  approximation  per  se  are  small 
as  well.  The  corresponding  ion-dynamical  corrections1  are  likely  to 
occur  only  near  h » « tu  from  line  center,  where  is  the  ion  plasma 

frequency.  For  the  plasmas  to  which  our  calculations  might  apply,  we 
have  iBp^  « 10  top,  and  by  comparing  this  with  the  average  Stark  shift 
from  Eq.  (h),  we  find  that  the  quasistatic  approximation  is  reasonably 
valid  for  all  densities  fulfilling 


a criterion  which  is  therefore  well  met  for  all  densities  exceeding  the 
low  - density  limits  for  the  utility  of  our  calculations  as  estimated 
by  Eqs.  (6)  and  (8). 

A corresponding  criterion  for  the  electrons,  obtained  by  using  the 
electron  plasma  frequency,  might  suggest  that  the  quasistatic  approxima- 
tion is  reasonably  valid  for  them  as  well,  in  spite  of  an  additional 
factor  10  in  the  condition  on  the  density.  However,  most  contributions 
to  the  field  are  characterised  by  frequencies  higher  than  tcp,  and  a 
more  quantitative  criterion  is  called  for  or,  conversely,  an  examination 
of  the  high  density  limit  for  the  validity  of  _ie  impact  approximation 
used  in  this  work.  This  limit  must  be  imposed  to  meet  the  requirement 
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y < 1 discussed  below  Eq,  (17) • Since  we  are  not  Interested  In  frequencies 
which  exceed  the  Inherent  level  splittings,  Eqs.  (19)  and  (17)  can  be 
used  to  find  that  y-values  of  Interest  fulfill 

y < z^/n^cT  < 1, 

unity  actually  being  approached  only  for  the  Ly-a  lines  and  near  the 
Inglls-Teller  limit.  (Note  that  nip  < & Is  always  true  near  the 
Inglls-Teller  limit  so  that  Is  Indeed  the  characteristic  frequency 
In  this  regime.) 

It  Is  clear  from  Eq.  (16),  that  the  collision  operator  becomes 

sensitive  to  the  actual  value  of  y as  this  parameter  approaches  unity. 

We  therefore  restricted  our  calculations  by  Imposing  the  requirement 

y £ 0.1,  In  which  case  a factor  2 uncertainty  In  y corresponds  to 

~ ± 30f>  error  In  the  approximate  collision  operator  while  the  excluded 

fraction  of  electrons  Is  entirely  negligible.  This  restriction  alone 

prevents  our  reaching  the  Inglls-Teller  limit  by  factors  < 10,  a gap 

that  could  be  narrowed  by  a (much  more  Involved)  unified  theory 
5 Ik  lc 

calculation.  * * ' [See  also  the  remarks  following  Eqs.  (2k)  and  (25).] 

Comparison  of  Eqs.  (12)  and  (16)  allows  the  identification 


^2S2  w y-1/2  M _25S 

pmin  ^min 


(20) 


w 

using j e Xdx/x  « jfcn(l/y)  and  replacing  impact  parameters  by  (dimenslon- 

y 

less)  angular  momenta,  L * mpv/fi.  According  to  Eq.  (13),  we  have 
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Given  the  restriction  y < 0.1,  we  therefore  need  to  consider  only  cases 

2 

with  L ^ 1.5  n and  should  expect  the  largest  total  error  if  L 

max  maty 


- is  near  this  limit.  (Errors  connected  with  p . will  be  almost  eliminated 

min 

by  the  modification  in  the  electron  broadening  operator  discussed  next) . 

IV.  ELECTRON  BROADENING 

. 

Since  the  diagonal  matrix  elements  of  the  collision  operator 
equal  one  half  to  the  total  rate  coefficients  for  electron  ion  (non- 
Coulomb)  collisions,  the  approximate  collision  operator  in  Section  II 
corresponds  to  total  cross  sections*’ 

*te)-  «> 

if  the  ion  is  in  the  n,£  level.  This  semiclassical  cross  section  for  the 
2p  and  l+p  levels  of  0 VIII  and  AjJXIII  is  compared  in  Table  I with  dis- 
torted wave  (quantum-mechanical)  calculations,1^  with  L = 20  in  both 

max 

cases  and  Lm^n  according  to  Eq.  (12)  for  the  semiclassical  cross  section. 
The  last  column  contains  a modified  semiclassical  cross  section,  in  which 
a strong-collision  term1'^  is  added  to  the  logarithm.  [fn(...)  is  re- 
placed by  + £n(...),  = 1.5,  = 0.75«]  After  this  modification 

the  distorted-wave  result  is  represented  by  Eq.  (22)  to  within  ~ 1C$, 
except  for  the  high  energy  values  for  2p,  The  distorted  wave  calculations 
as  such  should  be  almost  exact,  since  all  partial  wave  contributions  to 
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the  cross  section  stay  well  below  the  unitarity  limit  and  since  exchange, 
which  is  neglected,  is  not  likely  to  be  important  for  total  cross  sections 


(Corresponding  calculations  for  n = 3 and  5 give  = 1.0  and  = 0.5.) 

To  account  for  strong  (close)  collisions,  the  Maxwell-averaged 
collision  operator  estimated  by  Eq.  (16)  can  therefore  be  replaced  by 


9 = - ^ (fa) 1/2 » {if  l-i  (cn  ♦ I \ *"x  ?) 


(23) 


and  at  least  the  diagonal  matrix  elements  for  the  p states  may,  for 

y < 0.1,  be  expected  to  have  an  accuracy  of  ~ 20)6 . Comparisons  with 

IT  lft 

close-coupling  calculations  ’ for  2p  states  of  HI  and  He  II  support 
this  estimated  accuracy  for  the  Lyman-a  lines,  for  which  these  states 
are  most  important.  It  seems  reasonable  to  assume  that  Eq.  (23)  can  be 
used  for  all  levels  involved  in  our  calculations,  both  for  diagonal  and 
off-diagonal  matrix  elements  (in  the  parabolic  quantum  number  repre- 
sentation). While  we  suspect  that  the  strong  collision  terms  adopted 
here  may  lead  to  an  overestimate  for  these  off-diagonal  dipole- inter- 
action matrix  elements,  one  must  remember  that  the  corresponding  error 
in  the  calculated  profiles  may  be  counteracted  by  contributions  from 
off-diagonal  quadrupole- interaction,  etc.,  matrix  elements  which  were 
neglected . 

The  reader  may  also  wonder  whether  the  use  of  a Maxwell  distribution 
for  Eqs.  (16)  and  (23)  is  Indeed  appropriate  for  the  very  dense  plasmas 
in  which  Lyman-a  lines,  e.g.,  of  S1XIV  show  significant  Stark  broadening. 
To  address  this  question,  we  estimate  the  Fermi  energy  for  densities 
corresponding  to  the  Inglls-Teller  limit  in  Eq.  (10)  and  obtain 
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**=  £ o»v/5< 


(2k) 


This  characteristic  energy  is  therefore  considerably  smaller  than  typical 
thermal  energies,  kT  « z^E^A,  except  for  the  Lyman-a  lines,  for  which 
the  restriction  on  the  parameter  y prevents  us  from  calculating  profiles 
for  densities  close  to  this  limit.  Extensions  of  our  calculations  based 
on  the  unified  theory  to  these  densities,  should,  however,  account  for 
the  degeneracy  of  the  electron  gas  and  be  based  on  a fully  quantum- 
mechanical  treatment  of  all  multipole  interactions.  A corresponding 

19 

calculation  for  the  HI  Lyiuin-a  line  wings  now  exists  but  would  have  to 
be  generalized  to  allow  for  inelastic  collisions,  quasistatic  splitting 
and  degeneracy.  We  also  note  that  an  approximate  version  of  the 

unified  theory  suggests  corrections  to  the  impact  approximation  of  the 

-3  -1  21 

order  10  z for  y < 0.1.  An  application  of  the  unified  theory  to 

ion  lines  gives  corrections  of  order  lo’^nVzCz-l)2,  both  smaller  than 

those  corresponding  to  our  choice  of  the  maximum  impact  parameter.  (Note 

that  these  semiclassical  calculations  * are  based  on  the  dipole 

approximation.  They  would  have  to  be  supplemented  by  an  impact  parameter 

cutoff  near  the  excited  state  Bohr  radius  and  by  a quantum  calculation 

for  smell  L.) 


Lastly,  there  is  the  question  of  the  accuracy  of  the  semiclassical 
cross  section  of  ions  in  p states  for  partial  waves  L > 20.  We  there- 
fore compare  in  Table  II  the  L = 20  partial  wave  contribution  to  the 
cross  section  corresponding  to  Eq.  (22),  namely 


Lont  = 6"*i 


* (3)  ft) 


2.22...  1 

(n  -l  -f-1)  - , 


with  our  distorted  wave  quantum-mechanical  calculations.  Except  for  the 

lip  values  at  intermediate  energies,  the  agreement  Is  well  within  the 

20^  overall  theoretical  error  expected  for  our  calculations,  which  are 

presented  in  the  following  section.  (For  L 40,  deviations  were  found 

to  be  For  n = 2 and  3,  the  requirement  L p 20  imposes  according 

to  Eqs.  (20)  and  (21)  an  additional  restriction  (y  < 0.01  and  y < 0.05) 

on  the  parameter  y to  ensure  L values  large  enough  for  a quasiclassi- 

max 

cal  treatment  of  Debye  shielding,  finite  duration  of  collisions  and 
quasistatic  level  splittings.  For  the  Lyman-a  lines  we  must  therefore 
keep  the  density  below  the  Inglis-Teller  limit  by  factors  <50  or  treat 
inelastic  collisions  in  more  detail.  (Inelastic  collisions  are  respon- 
sible for  the  rather  large  deviations  between  semiclassical  partial 
cross  sections  and  distorted  wave  results  for,  say,  L = 10.) 


V.  RESULTS  AND  DISCUSSION 

Using  Eqs.  (11),  (17),  and  (23)  line  profiles  L(Aib)  were  calculated 
for  electron  temperatures  corresponding  to  near  maximum  abudnance  of  the 
ion  in  question  and  also  for  half  and  twice  these  temperatures.  Electron 
densities  were  varied  by  factors  of  ten  over  ranges  determined  as  dis- 
cussed in  Sections  I,  III,  and  IV.  Instead  of  the  angular  frequency  dis- 
placement from  line  center,  we  used  the  reduced  wavelength  separation 


0=  m = JMi.  10-a 


2»CF, 


with  the  Holtsmark  normal  field  strength  given  by  Eq.  (2). 
•8 

factor  10  is  connected  with  our  use  of  Angstrom  units.) 


(26) 


(The 

This  choice 


ensures  that  the  reduced  profiles 


If 


s(a)  = l(acd) 


d&a) 

da 


(27) 


have  a relatively  weak  dependence  on  electron' density.  From  Eqs.  (k)t 
(26),  and  the  Bohr  (Balmer)  formula  for  x follow  typical  o-values 

V 

1 p p _k  O 

a.  « 1.7  x 10  fg  n~  w 5 x 10  V (28) 

S 5 5 

e s z' 

for  a given  line  (with  principal  quantum  number  n of  the^ upper  state) 
and  radiating  ion  (of  nuclear  charge  z). 

Profiles  of  the  lines  of  Aj l XIII  are  Included  here  since  those 

22 

reported  previously  are  in  error  due  to  the  inadvertent  inclusion  of 
a section  of  code  that  adds  to  the  0 matrix  elements  (Eq.  12)  terms  to 
account  (approximately)  for faelastic  collisions.  In  conjunction  with  the 
Cq  of  equation  23  (equation  9 of  reference  22)  this  is  tantamount  to 
doubling  the  effect  of  inelastic  collisions. 

From  the  calculated  S(a)  profiles,  which  are  shown  in  Table  IV  and 
selected  lines  are  plotted  in  Figures  1-12,  we  determined  Stark  half, 
quarter,  and  eighth  intensity  widths  defined  by 


2 S(amax>  = S (al/2>* 

l S(°max>  = S <«l/4>’ 

5 S(0max)  = S (al/8)» 


(29a) 

(29b) 

(29c) 


where  a corresponds  to  the  maximum  of  S(a).  For  even  n-velues, 
max 

S(a)  has  a single  central  peak  (a  = 0),  for  odd  n-values  a central 

max 


minimum  and  two  symmetrically  located  maxima  (amaT ^ 0) « Table  III  con- 
tains values  of  a^,  Gfe,  o a^,  and  a^g. 

Also  given  in  Table  III  is  a quantity  ac  calculated  from 

°C  ” 5 (?  *)1/2°S  • «°> 

where  a is  the  parameter  defined  by  Eq.  (19)  and  ag  the  estimated 
Stark  broadening  from  Eq.  (28).  The  factor  1/5  corresponds  to  the 

estimated  ratio  between  Holtsmark  normal  field  strength  and  effective 

-•*  1 1/2 
mean  field  strength  used  in  Eq.  (5)  and  the  factor  (—  a)  is  from 

Eq.  (13)  of  ref.  11,  generalized  by  replacing  the  dipole  matrix  element 

2 

for  Lyman-a  (3  in  atomic  units)  with  n . 

According  to  ref.  11,  ac  is  a measure  of  the  broadening  from  field 
fluctuations  caused  by  electrons  in  the  Debye  shielding  clouds  of  ions. 

In  the  case  of  the  hydrogen  Lyman-a  line,  a substantial  theory-experi- 
ment10 discrepancy  (factor  ~ 2 in  half  width)  was  removed11  by  assuming 
a statistically  independent  Gaussina  broadening  with  1/e  width  given  by 
Q^,.  We  therefore  suggest  compounding  ac  with  the  thermal  Doppler  width 
before  convolving  our  calculated  Stark  profiles  with  Gaussian  profiles 
The  corresponding  modifications  are  only  important  near  a = 0,  but  do 
lead  to  substantial  modifications  of  central  profile  structures,  e.g.,  the 
narrow  peaks  of  a-  and  v-lines  and  the  central  dips  of  3-  and  6-lines. 

For  density  measurements,  reliance  on  these  central  structures  should 
therefore  be  avoided  until  the  physical  model  proposed  in  Ref.  11  has 
been  verified.  We  note  here  that  alternative  explanations  for  the  HI- 

Lyman  a discrepancy10  have  been  given  in  terms  of  ion-dynamical 

2k 

corrections  * which  would  probably  scale  differently. 
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Our  calculations  were  all  done  using  mlcrofleld  distributions 
calculated  for  the  case  of  equal  charges  on  radiating  and  perturbing 
ions,  i.e.,  for  p = z - 1.  As  can  be  seen  from  new  calculations  by 
Tlghe  and  Hooper,  ^ this  results  for  fixed  electron  density  in  profiles 
that  are  broader  than  those  one  would  obtain  if  most  of  the  perturbing 
ions  were  of  lower  charge  than  the  radiating  ions.  Also,  we  assumed 
equal  ion  and  electron  temperature  and  again  refer  the  reader  to  Ref.  23 
for  a discussion  of  distribution  functions  for  cases  where  these  temper- 
atures differ. 

On  the  other  hand,  although  our  calculations  are  strictly  valid 
only  for  one-electron  radiating  ions,  they  will  often  provide  a good 
approximation  for  resonance  lines  of  helium-like  ions  of  the  element 
with  the  next  higher  nuclear  charge.  This  will  be  true  as  soon  as 
densities  are  high  enough  for  the  profiles  to  be  broader  than  the 
electrostatic  splittings  of  np  and  ns,  etc.,  levels  of  the  unperturbed 
ions. 

Returning  to  Table  III,  it  is  clear  that  the  fractional  intensity 
widths  a1/n  for  a given  line  rarely  vary  by  more  than  a factor  of  two  1 
over  the  density  and  temperature  range  covered  by  these  calculations. 
Comparing  Ojyn  widths  for  analogous  lines  from  different  elements,  one 
finds  that  they  scale  approximately  as  z“5  for  near  maximum  abundance 
temperatures  and  densities  corresponding  to  similar  fractions  of  the 
largest  density  for  which  calculations  were  made.  This  scaling  is  in 
accordance  with  the  original  Holtsmark  theory^  and  therefore  with 
Eq.  (28).  It  suggests  that  reductions  of  the  quasistatic  broadening  by 
ions  due  to  ion- ion  correlations  and  Debye  screening  by  electrons  are 


largely  balanced  by  electron  Impact  broadening.  That  electron  impact 
broadening  remaina  sufficiently  Important  can  be  inferred  from  Eqs.  (k) 
and  (23).  which  with  kT  ~ z2  and  R»R  ~ z”2  indicate  that  the  ratio  of 
electron  and  ion  broadening  is  nearly  constant  along  the  iso-electronic 
sequence  if  the  electron  density  increases  as  z^.  This  tends  to  be  the 

case  under  typical  experimental  conditions  for  Stark  broadening  experi- 

26  27  2 8 

manta  ' ’ on  hydrogenic  ion  lines  in  the  Lyman  series. 
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Table  II 


• a) 

L-20  partial  wave  contributions  to  the  cross  sections  for  electron 

scattering  on  nr  levels  of  OVIII  and  A1XIII.  Semiclassical  values 

are  given  in  the  last  row. 


e/eh 

OVIII 

2p(xl02) 

4p 

e/eh 

AiXIII 

2p(xl03) 

4p(xl0) 

10 

2.02 

1.13 

20 

7.6 

4.1 

20 

1.83 

1.24 

50 

6.9 

4.7 

40 

2.15 

0.99 

100 

7.9 

3.9 

1.88 

0.98 

7.1 

3.7 

a)  2 

In  units  of  ita^  and  multiplied  by  E/E^  and  the  factors  following 
the  np  designations. 


Table  III 


Calculated  fractional  widths  of  Stark  profiles  for  Lyman  lines  from 

l*-n  transitions  in  one-electron  ions  for  various  temperatures  and  electron 
densities.  The  fractional  widths  are  in  units  of  10  ® Xngstrom  per  cgs 
field  strength,  as  are  the  positions  of  the  off-center  maxima  of 
the  n*3,  5,  and  7 lines  and  the  quantities  ac>  a^,  and  a,.,  which  are  defined 
in  the  text . 
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a 


0.5 

20 

0.00 

3 .60 

10.94 

1 .0 

20 

0 .00 

3.00 

15.47 

1.0 

21 

0.00 

3.70 

3.  33 

2.0 

20 

0.00 

2 .50 

21.88 

2 .0 

21 

0 .00 

3.10 

4.71 

2.0 

22 

0.00 

3.70 

1.  02 

CARBON 

1-3 

cr 

S » 78 

.0 

0.5 

19 

19.31 

4.70 

42.  85 

0.5 

20 

15.62 

5.70 

9.23 

1 .0 

19 

20  .25 

4.00 

60 . 60 

1.0 

20 

18.94 

4.80 

13.  06 

1.0 

21 

15.37 

5.90 

2 .31 

2 .0 

19 

20  .62 

3.40 

85.70 

2.0 

20 

20.06 

4.10 

18.46 

2.0 

21 

18.75 

4.90 

3.93 

2.0 

22 

14.87 

6 .00 

0.86 

CARBON 

1-4 

OC 

s = 123, 

0.5 

19 

0 .00 

7 .60 

40.63 

0.5 

20 

0.  00 

9.20 

8.  75 

1.0 

19 

0.00 

6.4  0 

57.46 

l .0 

20 

0 .00 

7.70 

12.38 

1.0 

21 

0.  00 

9.40 

2.  67 

2.0 

19 

0.00 

5.40 

81  .26 

2.0 

20 

0.00 

6.50 

17.51 

2.0 

21 
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Table  IV 

Stark  profiles  of  the  Lyman  lines  for  1 «-  n transitions  of 
one-electron  ions  for  various  temperatures  and  densities. 

The  profiles  are  tabulated  as  functions  of 


a 


1-^  1.  ■ (Angstrom  per  cgs  field  strength). 


Also  tabulated  are  the  parameters  Og,  ac,  end  Yg,  which  ere 
defined  in  the  text. 


The  theoretical  error  estimate  of  20^  must  be  increased 
(particularly  for  the  wings)  for  those  profiles  having 
yg  > 0.1  (0.01,  0.05  for  Lyman  a and  Lyman-p  respectively). 
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Vought  Corp.  (formerly  LTV  Aerospace  Corporation) 

Michigan  Division 

38111  Van  Dyke  Road 

Sterling  Heights,  Michigan  48077 

Attn:  Tech.  Lib. 


Ford  Aerospace  and  Communications  Corp. 

Ford  and  Jairi>oree  Roads 

Newport  Beach,  California  92663 

Attn:  Tech.  Information  Section 


Ford  Aerospace  and  Conmunications  Corp. 

3939  Fabian  Way 

Palo  Alto,  California  94503 

Attn:  D.  McMorrow  MS  G30 
Attn:  Library 


General  Electric  Company 
Space  Division 
Valley  Forge  Space  Center 
P.0.  Bo*  8555 
Philadelphia,  Pa.  I9IOI 

Attn:  J.  Peden  VFSC,  4230M 


General  Electric  Company  - Tempo 
Center  for  Advanced  Studies 
816  State  Street 
P.0.  Drawer  QQ 

Santa  Barbara,  California  93102 
Attn:  DAS1AC 


Institute  for  Defense  Analyses 
4 00  Army-Navy  Drive 
Arlington,  Virginia  22202 


Attn:  TDA  Librarian  R.  Smith 


47.  AVCO  Research  and  Systems  Group 
201  Lowell  Street 
Wilmington,  Massachusetts  OI887 

Attn:  Research  Lib,  A830  Rm  7201 


48.  BOM  Corporation 

7915  Jones  Branch  Drive 
McLean,  Virginia  22101 

Attn:  Technical  Library 


49.  Boeing  Company 
P.0.  Box  3707 

Seattle,  Washington  98124 
Attn:  Aerospace  Library 

50.  Dikewood  Industries,  Inc. 

1009  Bradbury  Drive,  S.E. 
Albuquerque,  New  Mexico  87IO6 

Attn:  L.  Davis 


51.  EG&G  Washington  Analytical  Services  Center,  Inc, 
P.0.  Box  10218 

Albuquerque,  New  Mexico  87114 
Attn:  Technical  Library 


52.  University  of  California 

Lawrence  Livermore  Laboratory 
P.0,  Box  808 

Livermore,  California  94550 

Attn:  L-18 

Attn:  L-153 

Attn:  John  Nickolls  A Div  L-545  (Class  L-33) 

Attn:  Tech.  Information  Dept  L-3 


53.  Sandla  Laboratories 
P.0.  Box  58OO 

Albuquerque,  New  Mexico  87H5 

Attn:  DOC  CON  for  3141  Sandia  Rpt.  Coll. 
Attn:  DOC  CON  for  5240  Gerald  Yonas 


mQm 


■ 


I 


5^ • Samso/SK 

Post  Office  Box  92960 

Vorldway  Postal  Center 

Los  Angeles,  California  90009 

(Space  Conm.  Systems) 

Attn:  SKF  Peter  H.  Stadier 


55  • AF  Weapons  Laboratory,  AFSC 

Kir tl and  AFB,  New  Mexico  87II7 

Attn:  CA 
Attn:  ELC 
Attn:  NT 
Attn:  SUL 
Attn:  DTP 


56.  Headquarters  USAF/RD 

Washington,  D.C.  20330 

Attn:  RDQSM 


57.  Samso/Dy 

Post  Office  Box  92960 

Worldway  Postal  Center 

Los  Angeles,  California  90009 

Attn:  DYS  (Technology) 


58.  Samso/lN 

Post  Office  Box  92960 

Worldway  Postal  Center 

Los  Angeles,  California  90009 

Attn:  Ind  MaJ  Darryl  S.  Muskln 
(Intelligence) 

39*  Samso/MN 

Norton  AFB,  California  92409 

Attn:  MNNH 
(Minuteman) 
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60.  Naval  Research  Laboratory 
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